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j. . Abstract. We explore the possibility to explain the LSND [1] result in the context of extra- 

' dimensional theories. If sterile neutrinos take shortcuts through extra dimensions, this results in 

altered neutrino dispersion relations. Active-sterile neutrino oscillations thus are modified and new 
types of resonances occur. 
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In theories with large extra dimensions, the Standard Model particles are typically 
confined to the 3 + 1 -dimensional brane, which is embedded in an extra-dimensional 
bulk [2, 3]. Singlets under the gauge group such as gravitons or sterile neutrinos however 
are allowed to travel freely on the brane as well as in the bulk. While the active neutrinos 
are confined to the brane, sterile states can take shortcuts through the extra dimension 
and the active- sterile neutrino oscillations [4] generate new resonances. 
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There are different ways in which these bulk shortcuts can be realized, one being the case 
of asymmetrically- warped extra dimensions. Warp factors shrink the space dimensions 
x parallel to the brane but leave the time and bulk dimension t and u unaffected [5, 6] 

dx 1 = dt z - e~ 2ku dx 2 - du 2 . (1) 

When the sterile neutrinos take shortcuts through an extra dimension a resonance arises 
due to an additional phase difference 8(Ht) = tdH + Hdt. Thus, in these models there 
are two sources of phase difference, the standard one tdH = LAm 2 /2E, and a new one 
Ht (St/t) arising from temporal shortcuts through the bulk available to sterile neutrinos. 
The two phase differences may beat against one another to produce resonant oscillations 
phenomena. 

By introducing the shortcut parameter e = (f brane - fbuik^brane _ § t j t the e ff ect i V e 
Hamiltonian for the two state system (one active and one sterile neutrino) is given by 

Am 2 / -cos20 sin20 \ „ £ / -1 
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Heff_ 4E { sin20 cos20 J E 2\ 1 (2) 



One notices that the diagonal terms in the effective Hamiltonian cancel out for a reso- 
nance energy E^ es = \J ~^^^ -- The effective mixing angle becomes 

. sin 2 26 

sin 2 26 = j- ( 3 ) 

sin 2 26+ cos 2 20 [l - ■£-] 

It can be seen that there are three distinct energy domains. Below the resonance energy 
one recovers vacuum mixing, at the resonance energy the effective mixing becomes 
maximal (6 = Jt/4), while above the resonance energy oscillations are suppressed. 

Starting from the asymmetrically- warped metric in Eq. (1), one needs to first solve the 
geodesic equations to calculate the time it takes the sterile neutrinos to travel through 
the bulk. For the sterile neutrinos to leave the brane, they must have a nonzero initial 
velocity uq along the extra dimension. When this happens, the geodesies are parabolic 
and they oscillate about the brane. The distance between two consecutive points where 
a geodesic crosses the brane is proportional to uq. Translational invariance is maintained 
on the brane and the Minkowski metric on the brane assures that Lorentz invariance 
is maintained on the brane. Therefore, momentum components are conserved on the 
brane, and we cannot generate a nonzero uq on the brane except as an initial condition. 
The uncertainty principle applied to the u dimension allows for such a nonzero velocity. 
Momentum conservation in the w-direction is a non-issue, as translational invariance in 
the ^-direction is broken by the brane itself. For a baseline L, as measured on the brane, 
there are more geodesic paths which cross the brane at the position of the detector, and 
the shortcut parameter for each of those paths is given by 

£„(v) = 1 - (") arcsinh Q , (4) 

where the scaling variable v = kL/2 was introduced, and the subscript n refers to 
neutrinos which enter the detector upon intersecting the brane for the n th time. 

These different modes have to be accounted for when calculating the probability of 
oscillation. Each mode is weighted by the quantum mechanical weight e' Sn , with S n 
being the classical action for the free particle. 

While the initial momenta are mostly on the brane, the uncertainty principle requires 
a nonzero p u component as well. Thus we assume a normalized GauBian distribution for 
the momentum component along the extra dimension with a width o which is related to 
the thickness of the brane. The distribution can be written in terms of the mode number 
through its dependence on velocity component uq. In order to be able to select only 
the geodesies which cross the brane at baseline length L, the integral over the momenta 
needs to be approximated with the corresponding sum with a measure An. 

When counting for all neutrinos An equals one, which is an upper bound. When 
looking only at the neutrinos which cross through the detector, the value of An is found 
by varying the action about the classical extremum. We account only for deviations 
AS = \S — S c \\ smaller than or of the order of h from the classical action, because larger 
variations lead to rapid oscillations and the integration averages to zero. 



Putting all the pieces together, the probability of oscillation including the weights 
mentioned above is given by 
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(n 2 + v 2 ) 3 /2 



with sin26>„ and 8H n obtained by replacing E Res with the resonant energy of the mode n. 
This now depends both on the energy E and the baseline L. The factor /3 is the velocity 
of the sterile neutrinos. 

Particularly interesting is the "Near Zone", defined for values v/n <C 1. A detailed 
motivation can be found in [7]. After making the small v/n expansion, a new feature 
which emerges is that the resonance peaks are functions of the energy and of the 
baseline through the combination LE rather than the energy alone as in the MSW matter- 
resonance [8, 9]. The LE dependence of the resonances is a novel feature of our model. 

The plot in Fig. 1 shows the probability of oscillation as a function of the baseline 
for the same value of energy. One can see the consecutive peaks corresponding to 
consecutive modes n. In our case the term sin(LdH n /2) oscillates fast, and phase- 
averaging then sets (sin ^4^-) to zero and (sin 2 to |. 

On the plot in Fig. 2, the dependence on the product LE of the oscillation probability 
becomes obvious. The resonance peaks are distributed in hyperbolic patterns on the L 
vs. E plane, with the peak closest to the origin of the axes being the one for n = 1. The 
relative height of consecutive peaks depends on the thickness of the brane o. Higher 
LE resonances are suppressed, and active-sterile neutrino mixing is suppressed for LE 
above the resonant values. 

The resonance encountered in our model might explain the observed excess in the 
LSND data. Sterile neutrinos decouple from the active neutrinos for long-baseline ex- 
periments as well as for high energies. Thus, no active-sterile mixing is expected in 
atmospheric data, in MINOS [10] or CDHS [11]. All explanations proposed so-far for 
the LSND and MiniBooNE anomalies assumed baseline-independent oscillations and 
mixing. Our model relies on metric shortcuts, and it does not discriminate between par- 
ticles and antiparticles. Thus it will be difficult to accommodate the recent MiniBooNE 
claims that an excess of flavor changing events exists in the neutrino channel [12] but 
not in the antineutrino channel [13]. It might still be possible though to explain the 
MiniBooNE data by non-standard matter effects [14]. The failure of previous models to 
reconcile short baseline data such as LSND with longer baseline data might be construed 
as favoring the extra-dimensional shortcut scenario. Finally, the bulk shortcut scenario 
might even relieve some of the remaining tension between the LSND and KARMEN 
[15] experiments since LSND has almost twice the baseline of KARMEN. 
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FIGURE 1. Oscillation probability as a function of the experimental baseline, for a GauBian distribution 
for uq (red and green curves). The green curve presents the phase-averaged oscillation probability, and the 
sinusoidal blue curve presents the probability as given by the standard 4D vacuum formula for oscillations 
between sterile and active neutrinos. Parameter choices are sin 29 = 0.003, k = 5/(10 8 m), E = 15 MeV, 
Am 2 = 64 eV 2 , and a = 100 eV. The resulting value of (LE)^ es is 550 m MeV. For our choice of E, the 
resonance peaks are found at the multiples L = n{LE)^/E = 37n m, n = 1,2,3 with the principal 
resonance corresponding to n = 1 . 



FIGURE 2. Oscillation probability (vertical) in the L-E plane, with the same parameters as in Fig. 1. 
Units of L and E are m and eV, respectively. 
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